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ABSTRACT 

Using the BSBM varying-alpha theory, with dark matter dominated by magnetic en- 
ergy, and the spherical collapse model for cosmological structure formation, we have 
studied the effects of the dark-energy equation of state and the coupling of alpha to 
the matter fields on the space and time evolution of alpha. We have compared its evo- 
lution inside virialised overdensities with that in the cosmological background, using 
the standard (A = 0) CDM model of structure formation and the dark-energy mod- 
ification, wCDM . We find that, independently of the model of structure formation 
one considers, there is always a difference between the value of alpha in an overdensity 
and in the background. In a SCDM model, this difference is the same, independent of 
the virialisation redshift of the overdense region. In the case of a wCDM model, espe- 
cially at low redshifts, the difference depends on the time when virialisation occurs and 
the equation of state of the dark energy. At high redshifts, when the wCDM model 
becomes asymptotically equivalent to the SCDM one, the difference is constant. At 
low redshifts, when dark energy starts to dominate the cosmological expansion, the 
difference between alpha in a cluster and in the background grows. The inclusion of 
the effects of inhomogeneity leads naturally to no observable local time variations of 
alpha on Earth and in our Galaxy even though time variations can be significant on 
quasar scales. The inclusion of the effects of inhomogeneous cosmological evolution 
are necessary if terrestrial and solar-system bounds on the time variation of the fine 
structure 'constant' are to be correctly compared with extragalactic data. 
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1 INTRODUCTION 

Studies of the time and spatial variations of the fine struc- 
ture 'constant', a, are motivated primarily by recent ob- 
servations of small variations in relativistic atomic struc- 
ture in quasar absorption spectra ( Murofiv ct ah] l200ll 
l2003l: I Webb et alj|200ll . [19991) . Three data sets, contain- 
ing observations of the spectra of 128 quasars continue to 
suggest that the fine structure 'constant', was smaller at 
redshifts z = 0.2 — 4.2 than the current terrestrial value 
ao = 7.29735308 x 10"^, with Aa/a = ^a(z) - qqI/qo = 
-0.57±0.11 X 10"^ i Murphy. Flambaum fc W ebb 2003). If 
this interpretation of these observations proves correct then 
there are important consequences for our understanding of 
the forces of nature at low energies as well as for the question 
of the links between couplings in higher dimensions and the 
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three-dimensional shadows that we observe llBarrow|l2002l) . 
Any slow change in the scale of the extra-dimensions would 
be revealed by measurable changes in our four-dimensional 
"constants". Further observational tests m ay help to pro - 
vide independent tests of the quasar data jDarhnell200ll . 
Attempts to generalise the standard model by including 
scalar fields which carry the space-time variation of the fine 
structure 'constant' could have important connections with 
the dark energy that is currently accelerating the expan- 
sion of the uni verse jParkinson. Bassett. fc Barrow I l2003l : 
IWettericbl2003f) and also create potenti ally detectable viola- 
tions of the weak equivalence p r inciple (|D vali fc Zaldarriagp 



12002'; 



Olive fc Pospelovl l2002l : iDamour fc ^olvakovr^& 



iNordtvec 



L200i 



Several theories have been proposed to in- 
vestigate the implications of a varying fine struc- 
ture 'constant'. Some based on grand u nifica- 

^ theories fbangacker. Segre fc Strassl^ |2iM 

bent fc FairbairnI 120031: iFritzschI 120021: ICalmet fc FritzschI 
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2002allbt iMarciand Il984l: iBanks. Dine fc Doudaj 120021: 
Armendariz-PicorJ 12002), some with extra-dimensions 



i BrsLX et al.' '2003':^PaIma et al.' "2003'; ' Youm" '2002'; 



Dend 



2002 ) 



Ml,- I Harrow, ^.a.nrtvik Magiienrl 

2002allbt iBarrow. Magueiio fc SandvikI l2002bl: 



12003 : .Paccetti Correia . Sc hmi dt fc Tayartkiladze 
others in four dimensions jBetensteil 
[Barrow. Sandvik fc Magueii " 



Sandvik, Barrow fc Magueiio 20o 3 iBarro^!^^^^^ 
2002 ; Huang & Li_ 2002 ; Kosteleckv. Lehnert fc Perry 
200d) and even in two dimensional black holes~ ilVagena£ 
2003r . All these models need to satisfy the present con- 
straints on the variations in a. These constraints can be 
divided into two main groups: local and astro-cosmological. 
The local constraints derive from experiments in our 
local bound gravitational system: the Oklo natural re- 
actor (z = 0.14) wher e < 10 "^ ("Fuiii et all |2000l: 
iDamour fc DvsonI Il996l: iFuiiil I2OO3I: Ishlvakhtor, hoO^ . 
However, the use of an equilibrium neutron spectrum in 
this analyses has recently been criticised by iLamoreauxl 
l2003l) . More realistic modelling leads to a best fit of the 
data to non-zero shift Aa/a = 4.5 x 10~*. ; the intra 
solar-system decay rat e ^^^Re ^^^Os, (z = 0.45), 
wher e < 10"'^ JOlive et al.l 120021: iFuiii fc Iwanotd 
l2003l: lOlive et all l2003l) : and the stability of terrestrial 
atomic clocks (z = 0), where — < 4.2 x 10 
iMarion et al.l 120021) . Other limits arise from 
equivalence exp eri ment s jDv^^^^^darri^ga| 
Ohve fc Pospelovl l2002l: iDamour fc Polvakovl 
Nordtvedd 120021: iMagueiio. Barrow fc SandvikI 
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but the limits they provide are more model dependent. 
In addition to the the above-mentioned positive signal 
from quasar absorption spectra, the astro-cosmological 
constraints are: the cosmic microwave background radiation 
(z = 10^), where | ^ | < 10"^ jMartins et al.l lioOS : 
Avelinoet al] I2OO0I: iBattve. Crittenden fc Welle: " 



Hannest adi I1999I: fKaplinghat. Scherrer fc Tur^ 



T 2001; 
I I1999D 



10"* 



10^ 



where 



and Big Bang nu cleosynthesis (z . 
1 ^ 1 < 2 X 10 -^ iMartins et a,l.ll20oi Ia^uo et ahlEoOfl 
iBattve. C rittenden fc Wellei' "2001*). Oth er em ission-line 
studies (fBahcall, Stoinhardt, & D . Schlcgell l2003l) give weak 
bounds of ^ < -2 ± 1.2 x 10"'' which are comparable 
to the earlier limits of refs. jCowie fc Songaiha Il995l: 



ifc ■ 

llvanchik. Potekhin fc Varshalovichlll999l: lLevshakovlll99l) 

There is a potential discrepancy between local and 
astro-cosmological constraints; in particular, between the 
constraint of Aa/a < at redshift < z < 0.45, comin i 
from the /3-decay in meteoritic samples ~lloiive et alJl200 
|20o3), and the explicit variation in a of Aa/a « 10~® at 
0.2 < z < 4.2, coming fro m the low-redshift end of the 
:)uasar absorption spectra (I Mu rphv. Webb fc F larnbaumI 



quasar absorption spectra ||MiyBnx^_WebpfeJL^ 
l2003l: I Murphv et aPboOll. l2003l : I Webb et al.ll200lL Il999 

A successful theory of varying a needs to explain this 
difference. This is a challenge. Models which use a very 
light scalar field, to drive variations in a, need an ex- 
treme fine tuning in order to satisfy the phenomenolog- 
ical constraints coming from geochemical data (Oklo, f3- 
decay), the present equivalent principle tests, and the quasar 
absorption spectra, simultaneously (Damour 2OO3I) . The 
presence of the dark energy plays an important role in 
any potential reconciliation because in simple scalar theo- 
ries iSandvik. Barrow fc Magueiidl2002^ any variation of a 
turns off as the universe starts to accelerate. 



A solution to thi s problem was proposed by us in 
jMota fc Barrowll200i) . where numerical simulations were 
shown to display behaviour for the evolution of a inside an 
overdensity that differs from that in the background uni- 
verse, during the formation of non-linear large-scale struc- 
tures. All previous studies of the variation of 'constants' 
in cosmology had neglected the effects of inhomogeneity 
in the Universe and the fact that our local observational 
constraints on varying constants are made within non- 
expanding matter overdensities. We found that the fine 
structure 'constant' evolves differently inside virialised clus- 
ters compared to the background universe. Specifically, a 
becomes a effectively time-independent inside a gravitation- 
ally bound overdensity after its virialisation. The fact that 
local a values 'freeze in' at virialisation, means we would 
observe no time variations in a on Earth, or elsewhere in 
our Galaxy, even though time-variations in a might still be 
occurring on extragalactic scales at a detectable level. For 
a typical galaxy cluster, the value of a today will be the 
value of a at the virialisation time of the cluster. Hence, 
the local constraints on time variations in the fine structure 
'constant', can easily give a value that is 10 — 100 times 
smaller than is inferred on extra-galactic scales from quasar 
absorption spectra. 

The reason why the growth and gravitational bind- 
ing of matter inhomogeneities affects the evolution of a 
is simple. Any varying-a theory implies the existence of 
a scalar field carrying the variations in a. This field is 
coupled to some selection of the matter fields, depend- 
ing upon their interactions jParkinson. Bassett. fc Barrow I 
l2003l:lAnchordoaui fc Go ldberg 2003; Wottorich 2003). Due 
to this coupling, any inhomogeneities of the latter will then 
affect the evolution of a. This effect is not only important 
in the non-linoar regime of large-scale structure formation 
{Barrow & O'Toolc 1999). Even in the linear regime of cos- 
mological perturbations, when these are small, ^ grows and 
tracks 



Sp„ 



^ i}^'^ during the dust-dominated era on scales 
smaller than the Hubble radius iBarrow fc MotallioO^ . It 
is therefore important to study the cosmological evolution 
of the fine structure 'constant' taking into account a more 
realistic universe, where matter inhomogeneities grow and 
lead to the formation of bounded objects. 

The dependence of the fine structure constant on the 
density of the matter fields leads to spatial variations of a. 
Spatial variations imply the existence of 'fifth force' effects. 
The fifth force induces an anomalous acceleration which de- 
pends on the material composition of the test particle, and 
so violates the weak equivalence principle (WEP). This is a 
general feature of any varying-a theory but he magnitude 
of the WEP violations is model dependent. In the case of 
the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) model, 
such violations are within the current WEP experiments 
(Barrow, M agueiio fc S andvik 2002a ). In part icular, it was 
shown in llMagueiio. Barrow fc SandvikI |2002|) that spatial 
variations of a in the vicinity of compact massive objects 
may well be within an order of magnitude below the exist- 
ing experimental bounds (see appendix IbI. 

There are two main features which affect the evolution 
of the fine structure 'constant' in a universe where large 
scale structures are formed. The most obvious, is the cou- 
pling between the scalar field, which drives variations in a, 
and the matter fields. The second is the dependence of non- 
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linear models of structure formation on the equation of state 
of the universe, and in particular that of the dark energy 
l|Lahav et alJll99ll : IWane fc Steinhardtl llQoi ) . In this pa- 
per we will study the dependence of a on these two quan- 
tities in the BSBM varying-a model. In this theory, vari- 
ations in Q are driven by the electromagnetically-coupled 
matter fields and the effects of inhomogeneity are large. 
In the next section we briefly summarise the BSBM model 
l|Sandvik. Barrow fc Mae ueiio 20o3 l and the spherical col- 
lapse model for the development of non-linear cosmological 
inhomogeneities. In section three, we investigate the depen- 
dence of the fine structure 'constant' on the equation of state 
of the dark energy. In section four, we show how spatial vari- 
ations in a may occur due to possible spatial variations in 
the coupling of a to the matter fields. We summarise our 
conclusions in section five. 



2 THE BSBM THEORY AND THE 
SPHERICAL-INFALL MODEL 

2.1 The Background 

We will study space-time variations o f a in the BSBM theory 
jSandvik. Barrow fc Magueiid l20o3) . which assumes that 
the total action is given by: 

S = y d^sV— 3 (-Cg + Cmatter + + CemS'^'^'^ (1) 

In this varying-a theory, the quantities c and h are taken 
to be constant, while e varies as a function of a real scalar 
field tl:{x'^), with e = eoe^ , = ^dy,i}}d^ip, ci; is a coupling 
constant, and Cem — —jfiiyf^"; i-matter is the Lagrangian 
of the matter fields. The gravitational Lagrangian is as usual 
jCg = — i^^q R, with R the Ricci curvature scalar, and we 
have defined an auxiliary gauge potential by = and 
a new Maxwell field tensor by /^^ = eF^,^ = d^a,^ — dvCi^, 
so the covariant derivative takes the usual form, = 9^ 4- 
ieoa^. The fine structure 'constant' is then given by a = 
aoe^^ with ao the present value measured on Earth today. 

The background universe will be described by a fiat, 
homogeneous and isotropic Friedmann metric with expan- 
sion scale factor a{t). The universe contains pressure- free 
matter, of density pm oc and a dark-energy fluid with 
a constant equation of state parameter, = p^/p^, and 
an energy-density p^ oc a"'^'^"'"™'*'. In the case where dark 
energy is the cosmological constant A, = pA = A/(87rG) 
and — —1. Varying the total Lagrangian, we obtain the 
Friedmann equation (?i = c = 1): 

H' = ^{p^{l + \C\e-'^) +p^+p,), (2) 

where H = d/a is the Hubble rate, = ^ip'^ , and = 
jC-em/pm is the fractiou of the matter which carries electric 
or magnetic charges. The value (and si gn) of C for baryonic 
and dark matter has been disputed ([ Dy^li Z^ldarri^gal 
l2002l:IOlive et alJl2002l:ISandvik. Barrow fc Magueiidl2002l) . 
It is the difference between the percentage of mass in 
electrostatic and magnetostatic fo rms. As explained in 
JSandvik. Barrow fc Magueiid l20o3) . we can at most esti- 
mate this quantity for neutrons and protons, with ("n ~ Cp ~ 
10~*. We may expect that for baryonic matter ^ ~ 10~*, 
with composition-dependent variations of the same order. 



The value of ( for the dark matter, for all we know, could be 
anything between -1 and 1. Superconducting cosmic strings, 
or magnetic monopoles, display a negative ^, unlike more 
conventional dark matter. It is clear that the only way to 
obtain a cosmologically increasing a in BSBM is with C < 0, 
i.e with unusual dark matter, in which magnetic energy dom- 
inates over electrostatic energy. 

The term |(^| e~^* represents an average^ of the (al- 
ways positive) energy density contribution from the non- 
relativistic matter which interacts electromagnetically. 

The scalar-field evolution equation is 

V; + 3Ji'V' = -- e"'*Cp- (3) 

In this article, we consider that the scalar field respon- 
sible for the variations of a is coupled to dark matter and 
to baryons. Both species are included in the C^pm terms. 
The coupling to dark matter appears in general in any dila- 
ton type theory, like string theory, where a scalar field is 
coupled to all the terms in the Lagrangian, but not nec- 
essarily in the same way. In the case of BSBM models, 
the coupling is motivated by certain clas ses of dark-matter 
models and their supersymmetric versions lOlive fc Pospelovl 
t2002J. The coupling between a and dark-matter could per- 
mit dark-matter-photon interactions of which would be me- 
diated by the scalar field. Due to the "invisible" nature 
of dark matter, its coupling to a is tightly constrained. 
In the case of B SBM models several con straints were im- 
posed on C in llOlive fc Pospelovl l2002l) . Other motiva- 
tions and conseque nces of this couphng are investigated in 
iBoehm et alJ l)2nn2h . where dark-matter-photon interactions 
are constrained using th e CMBR anis o tropie s and the mat- 
ter power spectrum. In iBoehm et al.l f2002) the cross sec- 
tion associated to the dark-matter - photon interaction is 
constraint to be a-y-DM /moM < 10~''ath/(100 Ge\/) « 
6 X 10~ cm^ GeV~^ , where aj-oM is the dark-matter- 
photon cross section, tudm is the dark- matter particle mass 
and (Jth is the Thomson cross section. Since in our case oth ~ 
o? ~ a^e*^ , then an approximate constraint to the coupling 
between the scalar field and dark-mater can be set imposing 
that a^-DMcr^-j < a^-oM, where a^-oM ~ 
and a^.-, ~ e"'*'''. Roughly giving |C| < lO'^ao ~ 10"^ 
when assuming the "usual" thdm ~ 100 GeV. In this article 
we will use |^| ~ 10~* as a reference value. In reality, as 
we show in section 4, the cosmological evolution of the fine 
structure constant is independent of value of |("| due to a de- 
generacy between the initial condition for tp and the coupling 
to the matter fields. The independence of our results with 
respect to the value of ^/uj give us the freedom to make the 
coupling to dark matter as small as desired. For instance, we 
are free to make the cross section photon - scalar field - dark 
matter smaller enough to satisfy all the constraints coming 
from the "invisible" nature of dark-matter. In a similar way, 
we are free to constraint the value of (^/ui in order to avoid 
problems associated with the so called self-interacting dark 
matter, which could radiate and cool, due to its 'magnetic' 

^ Each species in the Friedmann equation may have associated 
with it a different value of f . However, since we are interested 
in the large-scale cosmological behaviour we have averaged all 
matter contributions into a mean effective 'cosmological' value of 
C, (see appendix IXl. 
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nature. All these features could be combined to constraint 
the value of (. However that would be a lengthy and a highly 
model-dependent analysis, and will be considered by the au- 
thors elsewhere. 



2.2 The Density Inhomogeneities 

In order to study the behaviour of the fine structure 'con- 
stant' i nside overdensities w e will use the spherical-infall 
model jPadmanabhanl Il995ll . This will describe how the 
field in the overdensity breaks away from the field in back- 
ground expansion. The overdense sphere behaves like a spa- 
tially closed sub-universe. The density perturbations need 
not to be uniform within the sphere: any spherically sym- 
metric perturbation will evolve within a given radius in the 
same way as a uniform sphere containing the same amount 
of mass. Similar results could be obtained by performing the 
analysis of the BSBM theory using a spherically symmetric 
Tolman-Bondi metric for the background universe with ac- 
count taken for the existence of the pressure contributed by 
the dark energy and the field. In what follows, therefore, 
density refers to mean density inside a given sphere. 

Consider a spherical perturbation with constant inter- 
nal density which, at an initial time, has an amplitude Si > 
and \Si\ 1. At early times the sphere expands along with 
the background. For a sufficiently large Si, gravity prevents 
the sphere from expanding. Three characteristic phases can 
then be identified. Turnaround: the sphere breaks away from 
the general expansion and reaches a maximum radius. Col- 
lapse: if only gravity is significant, the sphere will then col- 
lapse towards a central singularity where the densities of 
the matter fields would formally go to infinity. In practice, 
pressure and dissipative physics intervene well before this 
singularity is reached and convert the kinetic energy of col- 
lapse into random motions. Virialisation: dynamical equi- 
librium is reached and the system becomes stationary: there 
are no more time variations in the radius of the system, R, 
or in its energy components. This phase determines the fi- 
nal pattern of variations in a, which becomes a constant 
inside the virialised region, otherwise the system would be 
unable to virialise jMota fc Barrow! |2003() . Meanwhile, de- 
pending upon the equation of state of the dominant matter 
field, a can continue to change in the cosmological back- 
ground. This behaviour naturally creates a situation where 
time variation of a on large cosmological scales is accom- 
panied by unchanging behaviour locally within galaxies and 
our solar system. It would explain the apparent discrepancy 
between the results of the quasars absorption spectra ob- 
serva tions and the /3-d ecay rate deductions from meteorite 
data iOlive et alJl2002l) . 

The evolution of a spherical overdense patch of scale 
radius R{t) is given by the Friedmann acceleration equation: 



R 
R 



{Pci^ (1 + ICl e-''^^) + 4p^. + (1 + 3u-^Jp^,) , (4) 



where pcdm is the density of cold dark matter in the cluster, 
P4,^ is the energy density of the dark energy inside the cluster 
and p^^ = ^ipc^ where tpc represents the scalar field inside 
the overdensity. We have also used the equations of state 
PVc = Pj/'ci Pcdm. = and p^^ = w^^p^^. 

In the cluster, the evolutions of pcdm and p^^ are 



given by 

V^c + S — V'c = e 

R LO 



Cpc 



Pcdm 



P<t>a 



.R 



-3^ (1 + w^J P,^,, 



(5) 
(6) 
(7) 



We will evolve the spherical overdensity from high red- 
shift until its virialisation occurs. According to the virial 
theorem, equilibrium will be reached when T = T, 
is the average total kinetic energy, and U is the average to- 
tal potential energy in the sphere. Note that we obtain the 
usual T = condition, for any potential with a power-law 
form [U oc R"), which includes our case. It is useful to write 
the condition for virialisation to occur in terms of the po- 
tential energies associated the different components of the 
overdensity. The potential energy for a given component ' x' 
can be calculated from its general form in a spherical region 
(Landau & Lifshitz 1975): 



27r / ptot(j>xr dr, 



(8) 
(9) 



where ptot is the total energy density inside the sphere, (j>x is 
the gravitational potential due to the Px density component. 

In the case of a ACDM model, the potential energies 
inside the cluster are: 

(11) 



Ug 



(10) 



-ttGpaMR, 
5 



(12) 



where Ug is the potential energy associated with the uniform 
spherical overdensity, Ua is the potential associated with A, 
and C/^^ is the potential associated with tpc- M — Mcdm + 
M^^ is the cluster mass, with 



M,dm = YPcd™(l + lCle~''^^'"' 



-)R^ 



4-K 



-Pi,.R 



(13) 
(14) 

(15) 



The virial theorem will be satisfied when 

Tvir ~ —-^Ug + Ua — 2^7^^, 

where T^ir = iMv^i^ is the total kinetic energy at virialisa- 
tion and v^i^ is the mean-square velocity of the components 
of the cluster. 

Using the virial theorem 11511 and energy conservation 
at the turnaround and cluster virialisation times, we obtain 
an equilibrium condition only in terms of the potential en- 
ergies: 

^UGiZv)+2UA{Zv)-U^^{z.,) = UG{zta)+UAizta)+Ui,^{zta), (16) 

where z„ is the redshift of virialisation and zta is the red- 
shift at the turnaround of the over-density at its maximum 
radius, when R = Rmax and J? = 0. In the case where ( — 
and = = we reduce to the usual virialisation condition 
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for AC DM models with no variation of a jLahav et alll99ll: 
IWang fc Steinhardtlll998l) . The generahsation from the cos- 
mological constant case to a dark-energy fiuid with a general 
equation of state is straightforward. One just needs to 
note that the potential energy associated to the dark-energy 
fluid is 

o 

where M^p^ is mass associated of that fluid, which is given 
by 

2.3 Setting the Initial Conditions 

The behaviour of the fine structure 'constant' during the 
evolution of a cluster can now be obtained numerically by 
evolving the background Friedmann equations J^J and ^ 
together with the cluster evolution equations @"Gl until 
the virialisation condition holds. In order to satisfy the con- 
straints imposed by the observations, we need to set up ini- 
tial conditions for the evolution. Since the Earth is now in- 
side a virialised overdense region, the initial condition for tp 
is chosen so as to obtain our measured laboratory value of a 
at virialisation, ac{zv) = a„ = ao- But, since the redshift at 
which our cluster has virialised is uncertain, we will choose 
a representative example where virialisation occurs over the 
range < z,, < 5. This is just for illustrative purposes, since 
in reality, the initial condition for t/i needs to be fixed only 
once, for our Galaxy. Hence, a in other clusters will have 
a lower or high er value (with respect to ao) depending on 
their Zv values jMota fc Barrowlboogh . 

After we have set Zv for Earth, another constraint 
we need to satisfy is given by th e quasar observations 
i Murphv. Webb fc Flambauml20o3) . This means that when 
comparing the value of the fine structure 'constant' on 
Earth, at its virialisation, q:„ = ao, with the value of the fine 
structure 'constant' of another region at some given redshift 
in the range accessed by the quasar spectra, 3.5 > z > 0.5, 
we need to obtain Aa/a = {a{z) — av)/cy.v ~ —5.4 x 10"''. 
This raises the question as to the location of the clouds 
where the quasar absorption lines are formed: are they in 
a region which should be considered as part of the back- 
ground or in an overdensity with somewhat lower contrast 
than exists in our Galaxy? Unfortunately, this question can- 
not be answered because we do not know the density of the 
clouds, only the column density. Nevertheless, these clouds 
are much less dense than the solar system. Because of this, 
it is a very good approximation to assume that the clouds 
possess the background density. 

Thus, the initial conditions for tp are chosen so as 
to obtain our measured laboratory value of a at viriali- 
sation adzy) = qq and to match the latest observations 
\ Murphv. Webb fc Flambauml20o3) for background regions 
at 3.5 >\z- z„| > 0.5. 

Here we may also wonder about possible measurements 
of a in Lyman-alpha systems and whether similar considera- 
tions should be applied because of local variations in density 
compared to this in the solar system. Assuming that the den- 
sity of the Lyman-alpha systems are much smaller than the 
Earth, we would expect to find a difference in a with respect 



to the value measured on Earth. This difference could even 
be of the same order as the one found in the quasar spectra. 
However, this is dependent on the density of the clouds and 
redshift we are making the measurements. Unfortunately 
we do not know the local density, only the column density 
along the line of sight. In order to make some numerical pre- 
dictions, we would need to create a model populated with 
clouds of different density but there are too many variables 
to make a reliable estimate of the effects at present. In the fu- 
ture it may be possible, with the accumulation of very large 
archives of data, to exploit the known differences in column 
density in the different types of system where the absorption 
lines are forms to search for correlations with column den- 
sity. The damped Lyman-alpha and Lyman-limit systems in 
the quasar absorption-system data sets have column densi- 
ties of order or exceeding 2 x lO^^cm"^ and 2 x lO^^cm"^ 
respectively. Again, rigorous exploitation of any perceived 
systematic trend in alpha between systems of different col- 
umn density would require the construction of a HI density 
model in the vicinity of the absorption system. 

3 THE DEPENDENCE ON THE 

DARK-ENERGY EQUATION OF SATE 

Non-linear models of structure formation will present dif- 
ferent featur es depending on the equation of state o f 
the universe llLahav et al.lll99ll: IWang fc Steinhardtlll99^ . 

The main difference is the way the parameter Ac = 
pcdm{zc) / ph[zv) cvolves with the redshift. The evolution of 
Ac depends on the equation of state of the dark-energy 
component which dominates the expansion dynamics. For 
instance, in a ACDM model, the density contrast, Ac in- 
creases as the redshift decreases. At high redshifts, the den- 
sity contrast at virialisation becomes asymptotically con- 
stant in standard (A = 0) CDM, with Ac ~ 178 at collapse 
or Ac ~ 148 at virialisation. This behaviour is common to 
other dark-energy models of structure formation (wCDM) , 
where the major difference is in the magnitude and the rate 
of change of Ac at low redshifts. 

The local value of the fine structure 'constant' will be 
a function of the redshift and will be dependent on the den- 
sity of the region of the universe we measure it, accord- 
ing to whether i t is in the background or an overdensity 
jMota fc Barrow! |2003|') . The density contrast of the viri- 
alised clusters depends on the dark-energy equation of state 
parameter, . Hence, the evolution of a will be dependent 
on as well. 

What difference we would expect to see in a if we com- 
pared two bound systems, like two clusters of galaxies? And 
how does the difference depend on the cosmological model 
of structure formation? These questions can be answered by 
looking at the time and spatial variations of a at the time 
of virialisation. The space variations will be tracked using a 
'spatial' density contrast, 

5a ^ Qc - otb j-^^^ 
a ab 

which is computed at virialisation (where ac{z — Zy) — a„). 
We assume there are no changes after this time. 

Since the main dependence of a is on the density of the 
clusters and the redshift of virialisation, we will only study 
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Figure 1. Value of a inside clusters as a function o/log(l + 2„), 
the epoch of virialisation. Solid line corresponds to the standard 
(A = 0) CDM, dashed line is the ACDM, dashed-dotted corre- 
sponds to a wCDM model with = —0.8, dotted-line corre- 
sponds to a wCDM model with = —0.6. In all the models, 
the initial condition were set in order to have ac{z = 0) = ao 
and to satisfy Aa/a ft: —5.4 X 10~® at 3.5 > \z — z^ \ > 0.5. 



Figure 2. The value of a in the background universe at the 
epoch of cluster virialisation, log(l + Zv). The solid line cor- 
responds to the standard (A = 0) CDM, dashed line is the 
ACDM , dashed-dotted line corresponds to a wCDM model with 
= —0.8, dotted line corresponds to a wCDM model with 
Wif, = —0.6. In all models, the initial condition were set in order 
to have ac{z = 0) = ao <*"<^ to satisfy Aa/a fx —5.4 X 10~^ at 
3.5 >\z- Zy\ > 0.5. 



dark-energy models where is a constant. Any effect 
contributed by a time-varying equation of state should be 
negligible, since the important feature is the average equa- 
tion of state of the universe. This may not be the case for the 
models where the scalar field responsible for the variations in 
a is c oupled t o dark energv l|P arkinson^ Bass ett. fc Barrow I 
ICor.ela,nd. Nnnes Posnelov I HSJ 
lAnchordoaui fc Goldberll2003D . 

In order to have a qualitative behaviour of the evo- 
lution of the fine structure 'constant', at the virialisation 
of an overdensity, we will then compare the standard Cold 
Dark Matter model, SCDM, to the dark-energy Cold Dark 
Matter, wCDM, models. In particular, we will examine the 
representative cases of = -1,-0.8,-0.6. All models 
will be normalised to have a„ = gp at z = and to sat- 
isfy t he quasar observations j Murphv. Webb fc FlambaumI 
I2OO3), as discussed above. This normalisation, although un- 
realistic (Earth did not virialise today), give us some indi- 
cation of the dependence of the time and spatial evolution 
in a on the different models. In reality, this approximation 
will not affect the order of magnitude of the spatial and time 
variations in a for the cases of virialisation at low redshift 
jMota fc Barrowll2003l) . 

3.1 Time-shifts and the evolution of a 

The final value of a inside virialised overdensities and its 
evolution in the background is shown in Figures and ^ 
respectively. From these plots, a feature common to all the 
models stands out: the fine structure 'constant' in the back- 
ground regions has a lower value than inside the virialised 
overdensities. Also, its eventual local value depends on the 
redshift at which the overdensity virialises. 

As expected, the equation of state of the dark energy 
affects the evolution of a, both in the overdensities and the 
background. A major difference arises if we compare the 
SCDM and wCDM models. In a SCDM model, the fine 



structure 'constant' is always a growing function, both in 
the background and inside the overdensities, and the growth 
rate is almost constant. In a wCDM model, the evolution of 
a will depend strongly on whether one is inside a cluster or 
in the background. In the wCDM background, a], becomes 
constant (independent of the redshift) as the universe enters 
the phase of accelerated expansion. Inside the clusters, 
will always grow and its value now will depend on the red- 
shift at which virialisation occurred. The cumulative effect 
of this growth increases as we consider overdensities which 
virialise at increasingly lower redshift. 

These differences arise due to the dependence of the 
fine structure 'constant' on the equation of state of the 
universe and the density of the regions we are measur- 
ing it. In a SCDM model, we will always live in a dust- 
dominated era. The fine structure 'constant' will then be 
an ever-increasing logarithmic function of time, a oc ln(t ) 
iSandvik. Barrow fc MagueiidlioO^: iBarrow fc MotalbOOSTl . 
The growth rates of «(, and Oc will be constant, since Ac 
is independent of the redshift in a SCDM model. In a 
wCDM model, dark energy plays an important role at 
low redshifts. As we reach low redshifts, where dark energy 
dominates the universal expansion, Ub becomes a constan t 
JSandvik. Barrow fc Magueiidliooa iBarrow fc M ota'2002ll 
but A c continues to increase, as wili~a7" ilMoto fc BarrcyM 
l2003h . The growth of a„ becomes steeper as we go from 
a dark-energy fluid with = —0.6 to the A-like case of 
W0 = —1. The intermediate situation is where = —0.8, 
due to the dependence of Ac on w^. 

Similar conclusions can be drawn with respect to the 
'time shift' of the flne structure 'constant', (Aa/a), at viri- 
alisation, see Figures 13 and |1| 
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Figure 3. Variation of Aa/a inside clusters as a function of 
the epoch of virialisation, log(l + Zv). The solid line corresponds 
to the standard (A = 0) CDM, the dashed line is the AC DM, 
the dashed-dotted line corresponds to a wCDM model with = 
—0.8, the dotted line corresponds to a wCDM model with = 
—0.6. In all models, the initial conditions were set in order to 
have ac{z = 0) = qq and to satisfy Aa/a fa —5.4 X 10~^ at 
3.5 >\z- 2„1 > 0.5. 



Figure 5. Variation ofSa/a as a function o/ log(l + Zt,). Solid 
line corresponds to the standard (A = 0) CDM, dashed line is 
the ACDM , dashed-dotted line corresponds to a wCDM model 
with ui0 = —0.8, dotted line corresponds to a wCDM model with 
= —0.6. In all the models, the initial conditions were set in 
order to have ac(z = 0) = ao and to satisfy Aa/a —5.4 X 10~® 
at 3.5 '>\z — Zv\> 0.5. 




log(1+z„) 

Figure 4. Variation of Aa/a in the background at the epoch of 
cluster virialisation, log(l + 2;„). The solid line corresponds to the 
standard (A = 0) CDM, the dashed line is the ACDM model, 
the dashed-dotted line is = —0.8, the dotted line (w^ = —OS). 
In all models, the initial conditions were set in order to have 
ac{z = 0) = «o and to satisfy Aa/a ft: —5.4 X 10~® at 3.5 > 

\Z - Zy\ > 0.5. 



Figure 6. Variation of da/ a as a function as a function of Ac- 
The single-point corresponds to the standard (A = 0) CDM, 
the dashed line is the ACDM , the dash-dot line corresponds to 
a wCDM model with = —0.8, dotted-line corresponds to a 
wCDM model with = —0.6. In all the models, the initial 
conditions were set in order to have ac{z = 0) = ao and to 
satisfy Aa/a -5.4 X lO"'' at 3.5 > \z ~ Zv\ > 0.5. 



3.2 Spatial variations in a 

Spatial variations in a will be dramatically different when 
comparing the standard CDM model with the wCDM mod- 
els. In the SCDM model, the difference between the fine 
structure 'constant' in a virialised cluster (a«) and in the 
background (a^) will always be the same, Sa/a ^ 5.2 x 10~^, 
independently of the redshift at which we measure it, Figure 
Again, this is because, in a SCDM model, Ac is always 
a constant independent of the redshift at which virialisation 
occurs. The constancy of 5a/a is a signature of the SCDM 
structure-formation model, and it may even provide a means 
to rule out the SCDM model completely if, when compar- 



ing the value of Sa/a in two different clusters, we do not 
find the same value, independently of z„. A similar result 
is found for the case of a dark-energy structure formation 
model at high redshifts where 5a/a will be constant. This 
behaviour is expected since at high redshift any wCDM 
model is asymptotically equivalent to standard CDM. 

As expected, it is at low redshifts that the difference 
between wCDM and SCDM emerges. When comparing 
virialised regions at low redshifts, 5a /a will increase in a 
wCDM model as we approach z = 0. This is due to an in- 
crease of the density contrast of the virialised regions, Ac, 
and the approach to a constant value of ai,, see Figure |S| In 
general, the growth will be steeper for smaller values of w^, 
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Figure 7. Variaiion of logio{a/a) as a function of log(l + Zy). 
Solid line corresponds to the Standard (A = 0) CDM, dashed- 
line is the ACDM , dashed-dotted corresponds to a wCDM model 
with = —0.8, dotted line corresponds to a wCDM model with 
= —0.6. In all the models, the initial conditions were set in 
order to have ac(z = 0) = ao O'l^d to satisfy A.a/a — 5.4x 10~® 
at 3.5 >\z- z„| > 0.5. 

although there will be parameter degeneracies between the 
behaviour of different models, which ensure that there is no 
simple relation between Sa/a and the dark-energy equation 
of state. Note that, independently of the structure formation 
model we use, tic/Qc at virialisation is always a decreasing 
function of time, as shown in Figure Q 



4 THE DEPENDENCE ON THE COUPLING 
OF a- VARIATION TO THE MATTER 
FIELDS 

The evolution of a in the background and inside clusters 
depends mainly on the dominant equation of state of the 
universe and the sign of the coupling constant C,/ui, which 
is determined by the t heory and the dark matter's iden- 
tity. As was shown in (|Sandvik. Barrow fc Magueiiol 120021 
iBarrow fc Motall2002^ ■ ah will be nearly constant for an 
accelerated expansion and also during the radiation era 
far from the initial singularity (where the kinetic term, 
p^, can dominate). Slow evolution of a will occur dur- 
ing the dust-dominated epoch, where a increases logarith- 
mically in time for C, < Q. When C, is negative, a will be 
a slowly growing function of time but a will fall rapidly 
even during a curvature-dominat ed era) for C, positive 
Sandvik. Barrow fc Magueiidliooi) . A similar behaviour is 
found for the evolution of the fine structure 'constant' inside 
overdensities. Thus, we see that a slow change in a, cut off 
by the accelerated expansion at low redshift, that may be re- 
quired by the data, demands that < in the cosmological 
background. 

The sign of is determined by the physical character 
of the matter that carries electromagnetic coupling. If it is 
dominated by magnetic energy then < 0, if not then C, > Q. 
Baryons will usually have a positive C, (although Bekenstein 
has arg ued for negativ e baryonic C, in ref. |Bckcnstoin 20025, 
but see (lDamouil2003l) ). in particular 10 ^ for neutrons 
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Figure 8. Evolution of a in the background (dashed-line) and 
inside a cluster (Solid line) as a function of log(l -|- z). For 
(/ui = -10"*, -4 X 10-8,-7 X 10-** and -IQ-^ . The lower 
curves correspond to a lower value of C/w. The initial condi- 
tions were set in order to have ac{z = 0) = ao o-^d to satisfy 
Aa/a = —5.4 X 10~^ at 3 > \z — Zy\ > 1, in the case where 
f = -7 X 10-**. 



and protons. Dark matter may have negative values of ^, for 
instance superconducting cosmic strings have ( m —1. 

In the previous section, we have chosen the sign of 
to be negative so a is a slowly- growing function in time 
during the era of dust domination. This was done in or- 
der to match the latest observa tions which suggest that a 
had a smaller value in the past (I Murphv et alJl200ll. l2003l : 
I Webb et all l200ll Il999f) . This is a good approximation, 
since we have been studying the cosmological evolution of 
a during large-scale structure formation, when dark mat- 
ter dominates. However, we know that on sufficiently small 
scales the dark matter will become dominated by a baryonic 
contribution for which > 0. The transition in the dominant 
form of total density, from non-baryonic to baryonic as one 
goes from large to small scales requires a significant evolu- 
tion in the magnitude and sign of C,/io. This inhomogeneity 
will create distinctive behaviours in the evolution of the fine 
structure 'constant' and will be studied in more elsewhere. 
It is clear that a change in the sign oi C^/uj will lead to a 
completely different type of evolution for a, although the 
expected variations in the sign of (/lo will occur on scales 
much smaller than those to which we are applying the spher- 
ical collapse model here. Hence, we will only investigate the 
effects of changing the absolute value of the coupling, IC/'^I, 
for the evolution of the fine structure 'constant'. 

From Figure |H] it is clear that the rate of changes in 
at and Uc will be functions of the absolute value of C/^- 
Smaller values of \(^/uj\ lead to a slower variation of a. A 
similar behaviour is found for the time variations in Aa/a, 
see Figure]^ 

The faster variation in a and Aa/a for higher values 
of jC/i^l is also a common feature for 5a/a, see Figure [THl 
This is expected. A stronger coupling to the matter fields 
would naturally lead to a stronger dependence on the matter 
inhomogeneities, and in particular on their density contrast, 
A,. 

In reality, the dependence of a on the coupling (^/uj has 
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Table 1. Time and space variations in a obtained for the corresponding redshifts of virialisation, Zv, for (^/w = 
— 10—^. We have assumed a ACDM model. The indexes 'b' and 'c', stand for background and cluster respectively. The 
italic and bold entries correspond approximately to the level of the Oklo and /3-decay rate constraints, respectively. 
The quasar absorption spectra observations correspond to the values of Aq/qI;,. The initial conditions were set in 
order to have aa{z = 0) = ctQ. 
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Figure 9. Evolution of h.a/a in the background (dashed-line) 
and inside a cluster (Solid line) as a function of log(l + z). For 
C/ui = -10-8,-4 X 10-8,-7 X 10-8 and -IQ-^ . The lower 
curves correspond to a lower value of C/u!. The initial condi- 
tions were set in order to have ac{z = 0) = ao o.'n.d to satisfy 
Aa/a = —5.4 X 10-^ at 3 > \z — Zv\ > 1, in the case where 
( = -7 X 10-8. 



a degeneracy with respect to the initial condition chosen 
for if). This is clear from the scale invariance of equation 
ijHJ under linear shifts in the value of tp —> + const and 
rescaling of ^/u and t. It is always possible, to obtain the 
same evolution (and rate of change) of a and Aa/a for any 
other value of |C/i^|; see for example Tables [112131 and HI 
where we have tabulated the shifts Aa/a, 5a /a, and time 
change, a /a, obtained in the clusters and in the background 
for various numerical choices of (/uj and virialisation redshift 
Zv. The observed Oklo and /3-decay constraints on variations 
in a are highlighted in italic and boldface, respectively. 

In the plots of this section, we have chosen the value 
(^/oj — —7 X 10-8 to be the one which satisfies the current 
observations, but we could have used any value of \(/uj\ be- 
cause of the invariance under rescalings. However, once we 
set the initial condition for a given value of (^/lo, any devi- 
ation from that value leads to quite different future varia- 
tions in a. This can occur if there are regions of the universe 
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Figure 10. Evolution of 5a/a as a function Sp/p. For C,/lli = 
-10-8, -4 X 10-8, -7 X 10-8 and -10-^. The lower curves cor- 
respond to a lower value of (^/u). The initial conditions were set in 
order to have ciciz = 0) = qq and to satisfy Aa/a —5.4 X 10~^ 
at 3 > \z — Zv \ > 1, in the case where ^ = —7 x 10^8. 



where the dominant matter has a difi'erent nature, and and 
do possesses a different value (and even sign) of i^/cu to that 
in our solar system. The evolution of a in those regions 
may then be different from the one that led to the value of 
ac{z — z„) — ao on Earth. 

It is important to note that, due to the degeneracy be- 
tween the initial condition and the coupling to the mat- 
ter fields, there is no way to avoid evolving differences in 
a variations between the overdensities and the background. 
The difference will be of the same order of magnitude as 
the effects indicated by the recent quasar absorption-line 
data. This is not a coincidence and it is related to the fact 
that we have normalised a{z = Zv) = ao on Earth and 
Aa/a = 5.4 x IQ-'^ at 3.5 > \z — z^] > 0.5. So, any varying- 
a model that uses these normalisations will create a similar 
difference in a evolution between the overdensities and the 
background, independently of the coupling C/lo. 

We might ask: how model-independent are these re- 
sults? In this connection it is interesting to note that even 
with a zero coupling to the matter fields (which is unre- 
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Table 2. Time and space variations in a for the corresponding redshifts of virialisation, z^, with ^/w = —10-^ .We 
assumed a KCDM model. The indices 'b' and 'c', stand for background and cluster respectively. The italic and bold 
entries correspond approximately to the Oklo and /3-decay rate constraints, respectively. The initial conditions were 
set in order to have cxc{z = 0) = oq. 
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Table 3. Time and space variations in a for the corresponding redshifts of virialisation, Zi,, for (^/o; = —10—^ .We 
have assumed a ACDM model. The indices 'b' and 'c', stand for background and cluster respectively. The italic 
and bold entries correspond approximately to the Oklo and to the /3-decay rate constraints, respectively. The initial 
conditions were set in order to have ac{z = 0) = oq- 

alistic), = 0, there is no way to avoid difference arising 
between the evolution of a in the background and in the 
cluster overdensities. While the background expansion has a 
monotonically-increasing scale factor, a{t), the overdensities 
will have a scale radius, R{t), which will eventually collapse 
at a finite time. For instance, in the case where C = 0, equa- 
tions and @ can be automatically integrated to give: 

ip(xa~^ , oc 7?"^ (18) 

The difference between those two solutions clearly increases, 
especially after the turnaround of the overdensity, when that 
region starts to collapse. As the collapse proceeds, the big- 
ger will be the difference between the background and the 
overdensity. Variations in a between the background and 
the overdensities are therefore quite natural although they 
have always been ignored in studies of varying constants in 
cosmology. 

5 DISCUSSION OF THE RESULTS AND 
OBSERVATIONAL CONSTRAINTS 

The development of matter inhomogeneities in our universe 
affects the cosmological evolution of the fine stru cture 'con- 
stant' (iBarrow fc MotJ2003l : lMota fc Barrowl20 03h. There- 
fore, variations in a depend on nature of its coupling to 



the matter fields and the detailed large-scale structure for- 
mation model. Large-scale structure formation models de- 
pend in turn on the dark-energy equation of state. This de- 
pendence is particularly strong at low redshift s, when dark 
energy dominates the density of the universe iLahav et alJ 
llQgiilWang fc Steinhardtll998^ . Using the BSBM varying-a 
theory, and the simplest spherical collapse model, we have 
studied the effects of the dark-energy equation of state and 
the coupling to the matter fields on the evolution of the fine 
structure 'constant'. We have compared the evolution of a 
inside virialised overdensities, using the standard (A — 0) 
CDM model of structure formation and dark-energy modi- 
fication {wCDM). It was shown that, independently of the 
model of structure formation one considers, there is always 
a spatial contrast, 5a/ 01, between a in an overdensity and 
in the background. In a SCDM model, 5a/a is always a 
constant, independent of the virialisation redshift, see Table 
|5] In the case of a wCDM model, especially at low redshifts, 
the spatial contrast depends on the time when virialisation 
occurs and the equation of state of the dark energy. At high 
redshifts, when the wCDM model becomes asymptotically 
equivalent to the SCDM one, 5a/ a is a constant. At low 
redshifts, when dark energy starts to dominate, the differ- 
ence between a in a cluster and in the background grows. 
The growth rate is proportional to |w</)|, see Tables |S1|7| and 
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Table 4. Time and space variations in a for the redshifts of virialisation, Zv, for C,/uj = —10—^^. We have assumed 
a AC DM model. The indices 'V and 'c', stand for background and cluster respectively. The italic and bold entries 
correspond approximately to the Oklo and /3-decay rate constraints, respectively. The initial conditions were set in 
order to have ac{z = 0) = ao. 



|H| These differences in the behaviour of the fine structure 
'constant', its 'time shift density contrast' (Aa/a) and its 
'spatial density contrast' (5a/a) could help us to distinguish 
among different dark-energy models of structure formation 
at low redshifts. 

Variations in a also depend on the value and sign of the 
coupling, C,/uj, of the scalar field responsible for variations 
in Q, to the matter fields. A higher value of |C/'^|, leads to 
a stronger dependence on the density contrast of the mat- 
ter inhomogeneities. If the value or sign oi C,/uJ changes in 
space, then spatial inhomogeneities in a occur. This could 
happen if we take into account that on small enough scales, 
baryons will dominate the dark matter density. The sign 
and value oi C,/uj will change, and variations in a will evolve 
differently on different scales. If there are no variations in 
the sign and value of C,/uJ, then the only spatial variations 
in a are the ones resulting from the dependence of the fine 
structure 'constant' on the density contrast of the region in 
which one is measuring a. At first sight, one might conclude 
that the difference between and Qt, is only a consequence 
of the coupling of a to matter. In reality this is not so. It 
is always possible to obtain the same results for any value 
\C,/uj\ with a suitable choice of initial conditions, as can be 
seen from the results in Tables Q |21 0] and |3] 

The results of this paper offer a natural explanation for 
why any experiment carried ou t on Earth (Fuiii ot al 20^ 
IPrestaee. Tioelker fc M alcki 1995: Sortais ct al. 2001), or in 
our local solar system (lOlive et aL .2002.1 . gives constraints 
on possible time-variation in a that are much stronger than 
the magnitude of variation that is consistent with the quasar 
observations ( Murphy, Webb & Flambaum 20Q^ on extra- 
galactic scales. The value of the fine structure 'constant' on 
Earth, and most probably in our local cluster, differs from 
that in the background universe because of the different his- 
tories of these regions. It can be seen from the Tables [S] 
|S| |7| and |H| that inside a virialised overdensity we expect 
Aa/of ~ — 10~^ for < 1, while in the background we 
have Aa/a ~ — 10~®, independently of the structure forma- 
tion model used. The same conclusions arise independently 
of the absolute value of the coupling (/to, see Tables 
andU 

The dependence of a on the matter-field perturba- 
tions is much less important when one is studying effects 



of varying a on the early universe, for example on the 
last scattering o f the CMB or the course of primordial 
nucleosynthesis llMajtins_et_alJ |200j ; lAvelino et al.l I2OOOI: 
iBattve. Crittenden Welled T20^F "^In the linear regime 
of the cosmological perturbations, small perturbations in 
a will decay or become constant in the radiation era 
jBarrow fc Motal l2003ll . At the re dshift of last scatter - 
ing, z = 1100, it was found in llMota fc Barrow! l2003l) 
that Aa/a < 10~^. In the background, the fine struc- 
ture 'constant^j_wilI_be_a^^^ the radiation 
era llSandvik. Barrow fc Magueiid 1200^1 so long as the ki- 
netic term is negligible. A growth i n value of a will occu r 
only in the matter-dominated era (iBarrow fc Motal l2002ll . 
Hence, the early-universe constraints, coming from the CMB 
and primordial nucleosynthesis, are comparatively weak, 
Aa/a < 10~^, and are easily satisfied. 
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APPENDIX A: BSBM EQUATIONS 

From the BSBM action Q the equation of motion for -0 
comes 

□0 = -e-^'^C,^ (Al) 

The right-hand-side (RHS) of equation lAH represents a 
source term for tp, which includes all the matter fields which 
are coupled to it. These include not only relativistic matter 
(like photons), but as well as non-relativistic one that inter- 
act electromagnetically. It is clear that jCem, vanishes for a 
sea of pure radiation since L^rn. ~ — B^^ /2 = 0. The 
only significant contribution to a variation in tp comes from 
nearly pure electrostatic or magnetostatic energy associated 
to non-relativistic particles. In order to make quantitative 
predictions we then need to know how much of the non- 
relativistic matter contributes to the right-hand-side (RHS) 
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Zy 


^li, X 10-^ 


X 10-7 


X 10-24 


S|e X 10-22 


a 


0.00 


-5.203 


0.000 


0.959 


1.409 


5.20 


0.06 


-5.229 


-0.266 


1.047 


1.537 


5.20 


0.12 


-5.254 


-0.517 


1.1S7 


1.670 


5.20 


0.38 


-5.349 


-1.463 


1.553 


2.281 


5.20 


0.46 


-5.374 


-1.717 


1.690 


2.481 


5.20 


0.54 


-5.399 


-1.962 


1.830 


2.687 


5.20 


1.07 


-5.534 


-3.313 


2.859 


4.198 


5.20 


2.12 


-5.720 


-5.175 


5.273 


7.744 


5.20 


3.00 


-5.834 


-6.315 


7.668 


11.261 


5.20 


4.92 


-6.012 


-8.102 


13.824 


20.305 


5.20 



Table 5. (A = 0) CDM: Time and space variations in a for the corresponding redshifts of virialisation, Zv, in the 
(A = 0) standard Cold Dark Matter model. The indexes 'b' and 'c', stand for background and cluster respectively. 
The italic and bold entries correspond approximately to the Oklo and to the /3-decay rate constraint, respectively. 
The quasar absorption spectra observations correspond to the values of Aa/a|i,. The initial conditions were set for 
C/w = —2 X 10-*, in order to have adz = 0) = ao- 



Zy 


^li, X 10-6 

fx 1" 


X lO-'' 


X 10-24 


2|c X 10-22 


5a X 10-6 

a 


0.00 


-5.440 


0.000 


0.332 


1.036 


5.44 


0.08 


-5.453 


-1.130 


0.385 


1.121 


5.34 


0.17 


-5.466 


-2.077 


0.440 


1.207 


5.26 


0.25 


-5.478 


-2.886 


0.497 


1.293 


5.19 


0.40 


-5.502 


-4.208 


0.613 


1.469 


5.08 


0.50 


-5.517 


-4.914 


0.691 


1.586 


5.03 


1.08 


-5.599 


-7.828 


1.217 


2.354 


4.82 


2.02 


-5.710 


-10.372 


2.238 


3.831 


4.67 


3.02 


-5.806 


-11.980 


3.518 


5.684 


4.61 


4.59 


-5.924 


-13.626 


5.876 


9.107 


4.56 



Table 6. lu = —0.6 CDM: Time and space variations in a for the corresponding redshifts of virialisation, Zv, in 
the w = —0.6 CDM model . The indices '6' and 'c', stand for background and cluster respectively. The italic and 
bold entries correspond approximately to the Oklo and /3-decay rate constraints, respectively. The quasar absorption 
spectra observations correspond to the values of Aa/a|(,. The initial conditions were set for i^/u) = —2 x 10-*, in 
order to have ac{z = 0) = ao. 



Zy 


^li, X 10-6 


^|cXlO-7 


^\bX 10-24 


fIcX 10-22 


Oi 


0.00 


-5.412 


0.000 


0.340 


1.078 


5.41 


0.09 


-5.427 


-1.587 


0.401 


1.148 


5.27 


0.18 


-5.441 


-2.857 


0.462 


1.220 


5.16 


0.34 


-5.468 


-4.777 


0.586 


1.370 


4.99 


0.42 


-5.481 


-5.527 


0.649 


1.447 


4.93 


0.52 


-5.498 


-6.362 


0.731 


1.549 


4.86 


1.09 


-5.587 


-9.516 


1.260 


2.239 


4.64 


2.01 


-5.706 


-11.997 


2.247 


3.595 


4.51 


3.38 


-5.841 


-13.948 


4.023 


6.125 


4.45 


4.92 


-5.953 


-15.297 


6.362 


9.510 


4.42 



Table 7. w = —0.8 CDM: Time and space variations in a for the corresponding redshifts of virialisation, Zy, for 
the w = —0.8 CDM model. The indices 'b' and 'c', stand for background and cluster respectively. The italic and 
bold entries correspond approximately to the Oklo and /3-decay rate constraints, respectively. The quasar absorption 
spectra observations correspond to the values of Aa/a\f,. The initial conditions were set for (^/oj = —2 X 10—*, in 
order to have ac{z = 0) = ao- 
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X 10^'' 


X 10-7 


^li, X 10-24 


f |c X 10-22 


— X 10-'' 

a 


0.00 


-5.382 


0.000 


0.361 


1.068 


5.38 


0.09 


-5.398 


-1.618 


0.427 


1.128 


5.24 


0.18 


-5.414 


-2.889 


0.493 


1.193 


5.12 


0.34 


-5.443 


-4.769 


0.623 


1.329 


4.97 


0.41 


-5.458 


-5.486 


0.688 


1.401 


4.91 


0.51 


-5.476 


-6.272 


0.773 


1.498 


4.85 


1.07 


-5.573 


-9.141 


1.303 


2.161 


4.66 


2.14 


-5.723 


-11.660 


2.508 


3.822 


4.56 


3.31 


-5.844 


-13.132 


4.056 


6.045 


4.53 


4.83 


-5.962 


-14.418 


6.408 


9.469 


4.52 



Table 8. A CDM: Time and space variations in a for the corresponding redshifts of virialisation, Zy, in the ACDM 
model. The indices 'b' and 'c', stand for background and cluster respectively. The italic and bold entries correspond 
approximately to the Oklo and /3-decay rate constraints, respectively. The quasar absorption spectra observations 
correspond to the values of Aa/a\i,. The initial conditions were set for (^/uj = —2 X 10—**, in order to have ac{z = 
0) = aa. 



of equation lAH . This can be parametrised by the ratio 
C = iCem/pm, where pm is the energy density of the non- 
relativistic matter. This non-relativistic matter, which inter- 
act electromagnetically, contributes to the Friedman equa- 
tion ^ as pm\C\s~'^^ . 

For protons and neutrons, ( can be estimated from the 
electromagnetic corrections to the nucleon mass, 0.63MeV 
and — 0.13Me«, respectively iPvali fc Zaldarriaeal l2002f) . 
This correction contains the contribution, which is al- 

ways positive, and also term of the form jua'^, where 
is the quarks' current (Sandvik, Barrow & Maguciio 2003) • 
Hence we take a guidance value of ~ 10^* for protons and 
neutrons. 

Using the parameter the fractio n of electric and mag- 
netic energies may then be written as JSandvikll2002l) : 



pm 



Pm 



(A2) 



where and are the electric and magnetic energies 
respe ctively. Using lA2t in equation lAH we have iSandvikI 



(A3) 



Since we are interested in the cosmological evolution of q, 
instead of using both parameters (^^ and we will use 
throughout this articles, the cosmological parameter, (, de- 
fined as C = - C^, which in the hmit where > 
is positive, and when ^ is negative. Note that, the 
cosmological value of has to be weighted, not only by the 
electromagnetic-interacting baryon fraction, but also by the 
fraction of matter that is non-baryonic. Hence the value and 
sign of C, depends strongly on the nature of dark matter to 
which the field might be coupled. 



APPENDIX B: THE BSBM MODEL AND WEP 
VIOLATIONS 

In BSBM the test-particle Lagrangian may be split as Lt = 
-Cm + e~'^^£,em- Variation with respect to the metric leads 
to a similar split of the stress-energy tensor, producing an 
energy density of the form p((l — Ct) + Ct6~^''')i s-^d soamass 



of — ("t) -f- C,te~'^^), (assuming electric fields dominate). 
In order to preserve their rati os of Ct = Ce.ml P test particles 
may thus be represented by iMaeueiio. Barrow fc SandvikI 
l2002h 



drm((l-0)-Cte-"^)[-3M.i''i1*^^^=^(Bl) 



where over-dots are derivatives with respect to the proper 
time r. This leads to equations of motion: 



2Cte 



-d'^i: = 



(B2) 



(1 - CO - Cte-^* 
which in the non-relativistic limit (with Ct ^ 1) reduce to 

dV 

where (j) is the gravitational potential. Thus we predict an 
anomalous acceleration: 

a- — 



\ OJTV J 



(B3) 



(B4) 



Violations of the WEP occur because C* is substance depen- 
dent. For two test bodies with Ci and C2 the Eotvos param- 
eter is: 



2\ai - 02] CslCi - C2 



0.1 + a2 



(B5) 



This can be written more conveniently as the product o f the 
following 3 factors ijSandvik. Barrow fc Magueiioll2002h : 



ttCp 



(a 



(B6) 



where E denotes Earth. If we take Cn ~ Cp ~ ICp ~ C^l = 
0(10"*) then for typical substances the first factor is ~ 
IQ-®. In order to satisfy the e xisting WEP vio lations exper- 
imental bounds l) Wil]il200ll: lNordtvedtlf2002h . we need to 
"play around" with the two l ast terms of equation (|B6t. If 
we u se the value assumed in iSandvik. Barrow fc Magueiid 
l2002h C/^ ~ -10"*. Hence, we need C = 0{1) to produce 
rj = 0(10"^'^), just an order of magnitude below existing 
experimental bounds. If we instead use <^/lo « —10"* we 
would need C = 0(10"*) to satisfy the same experimental 
bounds. The choice in the value of (^/ui just depends on the 
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nature of dark-matter, but this is beyond the scope of this 
paper. Nevertheless, as shown in this paper, any choice of 
^/u! will give us the same cosmological behaviour for Aa/a 
due to a degeneracy with the initial conditions for i/;. Hence, 
the results presented in this paper are not affected at all by 
any experimental constraint imposed by WEP violations. 
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